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ABSTRACT 

We present optical observations and 2D hydrodynamic modeling of an isolated shocked ISM cloud. 
Ha images taken in 1992.6 and 2003.7 of a small optical emission cloud along the southwestern limb 
of the Cygnus Loop were used to measure positional displacements of ~ O'.'l yr _1 for surrounding 
Balmer dominated emission filaments and 0'.'025 — '.'055 yr _1 for internal cloud emission features. 
These measurements imply transverse velocities of ~ 250 km s" 1 and ~ 80 - 140 km s _1 for ambient 
ISM and internal cloud shocks respectively. A lack of observed turbulent gas stripping at the cloud- 
ISM boundary in the Ha images suggests that there is not an abrupt density change at the cloud- 
ISM boundary. Also, the complex shock structure visible within the cloud indicates that the cloud's 
internal density distribution is two phased: a smoothly varying background density which is populated 
by higher density clumps. 

Guided by the Ha images, we present model results for a shock interacting with a non-uniform 
ISM cloud. We find that this cloud can be well modeled by a smoothly varying power law core with 
a density contrast of ~ 4 times the ambient density, surrounded by a low density envelope with a 
Lorentzian profile. The lack of sharp density gradients in such a model inhibits the growth of Kelvin- 
Helmholtz instabilities, consistent with the cloud's appearance. Our model results also suggest that 
cloud clumps have densities ~ 10 times the ambient ISM density and account for ~ 30% of the total 
cloud volume. Moreover, the observed spacing of internal cloud shocks and model simulations indicate 
that the distance between clumps is ~ 4 clump radii. We conclude that this diffuse ISM cloud is best 
modeled by a smoothly varying, low density distribution coupled to higher density, moderately spaced 
internal clumps. 



Subject headings: ISM: individual (Cygnus Loop) — supernova remnants 
dynamics — shock waves — hydrodynamics 
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1. INTRODUCTION 

The interaction of shock waves with the interstellar 
medium (ISM) such as those associated with supernovae, 
stellar winds, bipolar flows, H II regions, or spiral density 
waves is a fundamental process in interstellar gas dynam- 
ics and is key to understanding the evolution and struc- 
ture of the ISM. The highly nonlinear interaction be- 
tween supernova generated shocks and interstellar clouds 
is often not suited to analytic approaches but requires 
a multidimensional hydrodynamics study of the shock- 
cloud problem using high resolution methods. 

A hydrody namical study of a sho c ked I SM cloud 
was made by Klein. McKee. fc Colellal £l994, hereafter 
KMC94), who found that the cloud may be destroyed 
by a series of instabilities associated with the post- 
shock flow of inter-cloud gas past t he cloud. Earlier 
wor k on this problem includes that of iWoodwardl l)1976t) 
and | Nh) tman. Falle. & Gaskelll (1 1 982f) . More recently, 
Poludnen ko. Frank. & Blackmanl i|2002l) studied the role 
that internal cloud structure plays in the destruction of 
the cloud. 

For investigating shocked ISM cloud physics, the in- 
teraction between a supernova (SN) shock and low den- 
sity diffuse ISM clouds is of particular interest. Super- 
nova remnants (SNRs) shape and enrich the chemical 
and dynamical structure of the ISM which, in turn, af- 
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feet the evolution of subsequent SNRs. The details of just 
how SN generated shock waves interact with interstellar 
clouds are not well understood. 

There are several limiting factors in attempting to com- 
pare model simulations to observed SNR shock cloud in- 
teractions. While models can be viewed edge on and 
rotated in two or three dimensions, shocked interstellar 
clouds are viewed only in projection, which leads to a 
complex appearing shock structure due to multiple and 
overlapping shocks. In addition, one observes only a sin- 
gle epoch, i.e., a 'snapshot', of the interaction. These 
factors make it difficult to understand and model the 
time dependent kinematics and detailed dynamical pro- 
cesses of the interaction. Also, unlike how they are of- 
ten modeled, real interstellar clouds are neither cylindri- 
cal or spherical in shape nor sharp edged, with interiors 
very likely non-uniform in density. Furthermore, many 
shocked interstellar clouds are dense enough so that ra- 
diative losses, which can alter the overall dynamics of the 
shock-cloud interaction , are important (Mel lema et alJ 
2002l lFragile et al.l2004 . Finally, the inclusion of an cm- 
bedded magnetic field can drastically alter the dynamics 
of the interaction. For instance, a strong, ordered mag- 
netic field can suppress dynamical instabi lity growth pre- 
dicted by fluid dynamic al simulations l)Mac Low et alJ 
Hflfll iFragTle et al Jl200fih . 

In looking for an 'ideal' shock-cloud interaction, one 
would like to avoid many of the aforementioned effects 
and the Cygnus Loop supernova remnant affords several 
distinct advantages. Because of the remnant's large an- 
gular size (2.8° x 3.5°), low foreground extinction [E[B— 
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V] = 0.08 mag; iParked Il967t iFesen. Blair, k Kirshnerl 
1982), and wide range of shock conditions, the Cygnus 
Loop is one of the better locations for studying the ISM 
shock physic s of middle-aged remnants. At a distance of 
550^g ° pc (|Blair et al J 120051) . it has a physical size of 
27 x 33.5 pc. Located 8.5° below the galactic plane, the 
remnant lies in a multi-phase medium containing large 
ISM clouds with a hydrogen density of n = 5 — 10 cm" 3 , 
surrounded by a low er density inter-cloud component of 
n « 0.1 - .2 cm" 3 llDeNoverlll975l) . 

Recently, iPatnaude et alJ l)2002j) studied a small, iso- 
lated cloud along the southwest limb of the Cygnus Loop 
which met many of the desired cloud properties for in- 
vestigating shock-cloud interactions. This cl oud is rela- 
tively small (~ 2' in radius; 0.32 pc at 550 pc: lBlair et alJ 
2005), exhibits a fairly uncomplicated, line-of-sight inter- 
nal structure, and lies isolated from other shocked ISM 
clouds. Moreover, the shock-cloud interaction is domi- 
nated by non-radiative, or 'Balmer-dominated' filaments, 
indicating that the cloud-shock dynamics is not signifi- 
cantly affected by post-shock radiative losses. 

Here we present a new analysis of this small shocked 
cloud. Proper motion measurements and inferred shock 
velocities of individual filaments in and surrounding the 
cloud are presented. These results were used to estimate 
the initial conditions for hydrodynamical model simula- 
tions of a shock interaction with an unmagnetized, lumpy 
cloud. In §2, these new observations are presented as well 
as the technique used to measure the filament proper 
motions. Model parameter estimates are then discussed 
in §3. Our hydrodynamical models are presented in §4, 
where proper motion and density estimates are imple- 
mented in the model initial conditions. Model results 
are presented in §5, and they are compared to the south- 
west cloud in §6 with our conclusions in §7. 

2. OBSERVATIONS 

Narrow passband Ha images of the southwest region 
of the Cygnus Loop were obtained on 7 July 1992 and 
29 August 2003 using the MDM 2.4 m Hiltner telescope. 
For the July 1992 images, four 600 s Ha filter (FWHM 
= 80 A) exposures were acquired with a Loral 2048 x 
2048 front side illuminated CCD yielding a spatial res- 
olution of 0'.'343 pixel" 1 . Details of the 1992 observa- 
tions and subsequent data reduction can be found in 
IPatnaude et all p¥)l. 

Two 1000 s Ha filter (FWHM = 15 A) exposures were 
taken in August 2003 with a SITe 2048 x 2048 back side 
illuminated CCD with a resolution of 0'.'275 pixel" 1 . Us- 
ing IRAF 3 , the data were bias-subtracted, flat-fielded, 
and cosmic-ray hits were removed. The resulting 2003 
epoch image is shown in Figure Globally, the cloud's 
morphology is nearl y identical to that see n in the 1992 
images (Figs. 2 & 3 IPatnaude et alJl2002l) . but close in- 
spection between the two epoch images showed measur- 
able proper motions for both internal cloud structures 
and the surrounding thin shock front filaments. 

3. ANALYSIS AND RESULTS 

3 IRAF is distributed by the National Optical Astronomy Obser- 
vatories (NOAO), which is operated by the Association of Univer- 
sities for Research in Astronomy, Inc. (AURA) under cooperative 
agreement with the National Science Foundation. 



Following the procedure de scribe d by 
Thorstcnscn. Fes en. & van den BerghJ |200l . the 
coordinate systems of the two Ha images were aligned 
using DAOPHOT in conjunction with the USNO-A2.0 
catalog. The datasets were then rebinned to an effective 
image scale of O'/l pix -1 . This rebinning introduced 
a small global offset of — 0'/07 pix -1 between the two 
images, uniform across the entire field of view. 

3.1. Proper Motion Measurements 

Individual filament regions for the proper motion anal- 
ysis were selected based on their projection onto the 
plane of the sky, the complexity of the filament and sur- 
rounding regions, and the brightness of the filament fea- 
ture. Based on these criteria, 14 filaments within the 
cloud, including both Balmer dominated and radiative 
filaments, and 21 regions from the surrounding Balmer 
dominated shock front were chosen (see Figure Q}. 

One-dimensional intensity profiles were extracted for 
each region and the pixel shift in each shock filament was 
computed usin g the IRAF task xcsao. which is based on 
the software of iTonrv fc Davia l|1979D . While this task 
was written to compute relative radial velocities via the 
cross-correlation function between two spectra, the cross- 
correlation function yields accurate filament motions in 
terms pixel shifts between two images. For thin Balmer- 
dominated filaments and bright and sharp cloud shock 
features, the cross-correlation analysis was able to match 
the shock fronts between the two epochs and measure the 
pixel shift between the two data sets to an accuracy of 
10% - 15%. The results from this analysis are listed in 
terms of proper motion (mas yr ) and transverse veloc- 
ity (km s" 1 ) in Tabled The quoted velocities assume a 
distance of 5501^° P c lIBlair et a,l Jl2"f)05l) . 

Example data and cross-correlation functions for two 
regions are shown in Figure [5] These one-dimensional 
filaments and cross-correlation functions are representa- 
tive of the data for the non-radiative filaments (Fig. 
left) and internal cloud filaments (Fig. |2J right) , where 
there is often a plateau of emission (from shocked cloud 
material) downstream from the shock front and then a 
steep rise in emission at the cloud shock front. 

Quoted errors in Table ^ include the signal to noise 
in the filament region, the curvature of the filament re- 
gion, the profile of the filament when convolved with the 
image PSF (FWHM 1992 = l'.'O; FWHM 20 o 3 = 0'.'7), and 
the contribution from the background, local nebulosity, 
and adjacent filaments. For well resolved filaments, the 
cumulative effect of these errors is ~ 0.6 pixel, or about 
10 km s -1 . 

3.2. Estimate of Cloud Parameters 

As discussed in IPatnaude et al.1 l)2002|) . this shock- 
cloud interaction is nearly tangent along the line of sight 
to the southwestern limb of the southern region of the 
Cygnus Loop. The cloud is being run over by the rem- 
nant's shock front which is moving roughly east-west. We 
have divided the cloud into two regions: the ISM shock 
front and the cloud-shock region. Based on the filament 
velocities listed in Tabled we estimated an interstellar 
shock velocity of 250 km s" 1 associated with the Balmer- 
dominated filaments. The wide range of Balmer filament 
velocities observed (140 — 260 km s" 1 ) may be due in part 
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to density fluctuations around the cloud and the fact that 
only one component of the filament velocity is measured. 
That is, for filaments which are highly curved, the space 
velocity of the shock might be 200 km s _1 , but the lo- 
cal measured velocity might be in a direction other than 
perpendicular to our line of sight. Furthermore, though 
filaments were chosen based upon selective criteria, fac- 
tors such as low signal to noise as well as adjacent, over- 
lapping filaments contributed in some cases to a poorer 
cross-correlation between the two images. Nonetheless, 
our estimated shock velocity ~250 km s _1 is consistent 
with the X-ray shock velocity of ~ 300 km s -1 inferred 
from the ROSAT PSPC measurements l|Patnaude et alJ 
2002). 

The shocked cloud can be further divided into regions 
where the cloud-shock is interacting with the cloud, and 
where it is interacting with cloud clumps or "cloudlets" . 
In general, the inferred shock velocities vary widely (65- 
140 km s _1 ). This suggests that the density structure of 
the cloud is fairly complex, as the cloud-shock appears to 
have been slowed less in certain areas relative to others. 

Based on these measurements, we adopt a shock ve- 
locity range inside the cloud of 60 — 100 km s _1 . These 
estimated cloud-shock velocities in turn imply a range of 
density contrasts in the cloud. Assuming ram pressure 
equilibrium, (p a Vg ~ PcV 2 s ) the density contrast between 
the cloud and the ambient medium, x = Pel Pa, is ~ 4- 
17, with higher values representing areas populated with 
cloudlets, and lower values representing regions of low 
density within the cloud. 

The low density nature of this cloud permits us to view 
its internal shock dynamics. We have used the structure 
and spacing of the internal shock fronts to estimate the 
dumpiness of the cloud. The easiest place to do this is at 
the western nose of the cloud (Regions C8-C10). Mea- 
surements suggest that the post-shock spacing of clumps 
in the cloud is ~ 10" - 30". The upper limits corresponds 
to shocks which are more highly evolved, while the lower 
limit corresponds to 'small' shocks. Based on the size of 
the cloud (2' - 4' diameter), the cloudlets likely account 
for 30% of the total cloud volume. Furthermore, based 
on a maximum compression of ~ 4 for cloud clumps, 
we estimate the spacing between cloudlets to be w 4-5 
cloudlet radii (a c i ou diet) ■ 

4. HYDRODYNAMICAL MODELS 

For modeling this shock-cloud interaction, we have 
used the numerical hydrodynamics soft ware VH-1 
ijBlondin fc Knerrlll992t Istevens et al.lll992|) . which im- 
plements the piecewise parabolic method (PPM) to solve 
the e quations of gas dynamics ijColella. fc Woodward 
1984). The PPM approach incorporates a fixed computa- 
tional grid to evolve the standard conservation equations 
of mass, momentum, and energy. We assume an ideal, 
inviscid fluid with a ratio of specific heats, 7, equal to 
5/3. VH-1 does not explicitly treat the collisionless shock 
physics associated with Balmer-dominatcd shocks. How- 
ever, the goal of this study is to understand how a blast 
wave interacts with a diffuse ISM cloud. For these pur- 
poses, VH-1 serves as an excellent tool for tracing the 
motion and dynamics of this interaction. 

Simulations were performed on a 1440 x 1440 Carte- 
sian grid. The fiducial physical size of the square grid is 
lxl, with dx = dy = 6.9 x 10~ 4 , The size of the cloud 



sets the scale of the models. The cloud is w 2' in radius 
(w 10 18 cm x dsso pc ). On average, the cloud radius is 
30% of the computational grid, or w 450 cells per cloud 
radius. Therefore, the physical length scale of the grid 
Ax w 2 x 10 15 cm cell" 1 . 

We estimate the importance of radiative losses by cal- 
culating the cooling time scale and comparing it to dy- 
namically relevant time scales (mainly, the cloud crush- 
ing time and the pressure variation time scale). In gen- 
eral, radiative losses will be considered important if the 
cooling time is comparable to or shorter than the cloud 
crushing time. We estimate t he cooling time using the 
approximation of Kahnj l|1976[) . t coo i — Cv 3 /p, where v 
is the cloud shock velocity in units of km s -1 , p is the 
cloud density in units of gm cm -3 , and C is a constant = 
6.0 x 10~ 35 . Assuming a cloud shock velocity of 140 km 
s _1 and a cloud density of 10~ 24 gm cm~ 3 , we estimate 
a cooling time t coo i > 5000 yr. In contrast, the cloud 
crushing time, t cc = X 1 a o/vbi is ~ 3800 yr, for a blast 
wave velocity of 250 km s _1 and an initial cloud radius 
of cto = 0.3 pc. Furthermore, the pressure variation time 
scale is ~ 0.1i cc (KMC94), which is < t coo z- Thus, it ap- 
pears that neglecting the effects of cooling in our models 
will not have a significant impact on our results. This 
is supported bv lPatnaude et all lj2002j) who showed that 
this cloud is only weakly radiative. 

Under the assumptions that radiative losses are not 
dynamically important and that magnetic fields are not 
present (or that the cloud is only weakly magnetized) , the 
shock-cloud interaction can be wholly defined by two pa- 
rameters (KMC94), the shock Mach number, and x, the 
density contrast of the cloud. Assuming an ISM sound 
speed of « 10-15 km s _1 and a blast wave velocity of 250 
km s _1 , we estimate a shock Mach number of M ~ 20. 
Furthermore, as pointed out in §3, we estimate a cloud 
to ISM density contrast of x ~ 6. 

To further simplify the problem, we chose a set of non- 
dimensional variables such that the ISM density p a is 
set to the ratio of specific heats, 7 = 5/3, and the ISM 
pressure, P a , is set to unity. The ISM sound speed, c a , is 
thus set to 1 (c a = (7 P a j Pa) 1 / 2 ), and the shock velocity 
Vb is just the shock Mach number. 

Model results by KMC94 suggested that a cloud radius 
should be at least of order 120 cells. For our models here, 
we chose the main cloud to have a radius of 300-500 cells 
(~ 10 17 ~ 18 cm). The internal cloudlets have sizes which 
are 10-20% of the cloud radius. Therefore, the cloudlets 
are only 33% the suggested size. This small cloudlet 
size limits our ability to resolve instability growth along 
their boundaries. However, the goal here is to understand 
the global, internal morphology of the cloud, rather than 
small scale mixing within the cloud. 

We broke the cloud's density structure into two parts: 
A background density profile, and a clumping or density 
perturbation distribution. The background density dis- 
tribution is chiefly responsible for the large-scale shock 
features, such as how the shock drapes over the cloud 
edges, while also defining the initial internal cloud shock 
velocity. In contrast, the internal cloud density pertur- 
bations lead to the formation of small scale shock struc- 
tures within the cloud and have little effect on the cloud- 
ambient medium boundary layer. 

Based on the cloud's emission features (Fig.^J, we as " 
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sume that the large scale density structure of the cloud 
is smoothly varying. The interface between the blast 
wave and the cloud shock, seen along the northern and 
southern periphery of the cloud, appears smooth. This 
suggests that the cloud is surrounded by a low density 
envelope. There is no evidence to suggest that the cen- 
tral density is sharply peaked, so we assume that at some 
inner radius the density profile turns over and becomes 
relatively constant throughout. Therefore, we assume 
that the cloud consists of a smoothly varying core sur- 
rounded by a low density envelope. A function which 
fits this description is a truncated Lorentzian coupled to 
a power law core: 



where r 2 = a^Kxmax — 1), A = 0.11 sets the maximum 
core density, and a is the power law index = 0.5 which 
ensures continuity across the core-envelope interface. 

While the fine-scale structure of density perturbations 
within the cloud is not known, the lack of observed 
dynamical instabilities (at the resolution of the obser- 
vations), suggests that such perturbations are proba- 
bly smoothly varying. Therefore, for models using the 
above cloud density distribution, we chose to model the 
cloudlets as Gaussians. The spacing of the Gaussians is 
such that Ar « 4er between the cloudlet cores, which is 
consistent wit h the optimal spacing of Ar ~ A.2 a ^. ml ,dip . t 
suggested by Poludnenko, Frank. & Blackman (2002). 
Individual model parameters are listed in Table |21 For 
comparison, we include models of cylindrical clouds with 
similar density distributions. 

5. RESULTS 

Our basic shock-cloud interaction, Model 1, is shown 
in the top panel of Figure [3] This model is of a Mach 
20 shock interacting with a cloud of radius ao — 0.3 and 
constant density contrast x — 6. Figure shows the 
model at t = 1.4, 2.6, 3.6, and 4.7 x 10~ 2 . Panel c 
(t = 3.6 x 10 -2 ), shows the model at approximately one 
cloud crushing time (the cloud crushing time, given by 
Equation 2.3 of KMC94, is t cc w 3.7 x 10~ 2 ). According 
to KMC94, the growth time for Kelvin-Helmholtz (K-H) 

instabilities is txH ~ X 1 ^ 2 /(kv re i), where v re i is the rela- 
tive velocity between the shocked cloud and the shocked 
ambient medium. 

The relative velocity between the post cloud-shock 
material and the post shock ambient material is given 
to first order by the relative jump conditions between 
the cloud and the ambient medium. Since M oc x 1 ^ 2 ! 
i>KH ~ X 1 ^ 2 tcc (for k ~ ao). Higher wavenumber pertur- 
bations will form on a shorter time-scale. In Panel c of 
Figure |3 there is clear evidence for K-H growth along 
the backside of the cloud. In fact, it is evident that K-H 
growth occurs much earlier (top, Panel 6, Fig. [SJ. 

Model 2 is shown in the bottom panel of FigureOl This 
model has the density distribution described in Equa- 
tion ^ with an inner radius r*j = 0.6ao. The evolution 
of this model is markedly different than that of Model 1 
(Fig El . The \ = 10 n sted in Table is what the den- 
sity would be at the center of the Lorentzian. However, 
since the Lorentzian is truncated at r,, the effective x 1S 



much lower, by an amount 1/(1 -fr 2 /^). Therefore, the 
maximum x i n the cloud is ~ 3, or half the x of Model 
1. More importantly than the lower x between the two 
models, Model 2 does not show signs of the instability 
growth seen in Model 1. This interesting result lends 
weight to the notion of a smooth boundary between the 
ISM and an embedded ISM cloud. 

While Models 1 and 2 accentuate the differences which 
arise between a smoothly varying density distribution 
and that of a sharp edged cylinder, the remaining models 
(3-6) simulate how the internal density structure affects 
the cloud shock. Models 3 and 4 (Fig. 0J are of cylindri- 
cal clouds of x = 6 and 3. Both clouds contain perturba- 
tions with a x of 10 above the cloud density (or 30 and 
60 times the ambient density). In Model 3 (x = 6), the 
flow around the cloud is still strongly influenced by the 
higher cloud density. The inclusion of cloudlets results 
in the formation of shock structure within the cloud (not 
seen in Model 1). However, in the late time (Panel d), 
the morphology of the cloud is still similar to the late 
time morphology of Model 1 . 

The evolution of Model 4, on the other hand, is more 
strongly influenced by the cloudlets. This is because 
the density contrast between the cloud and the ambient 
medium here is only 3, and thus the cloud shock velocity 
is not significantly different than the blast wave velocity 
(V3 lower). Instead, the blast wave is more influenced 
by the high density cloudlets (relative to the cloud), as 
seen in Figure 0] (bottom). While this model reproduces 
the observed shock diffraction (c.f. Fig. the presence 
of instability growth along the (albeit low density) model 
cloud boundary is not something observed in the obser- 
vations. 

Models 5 and 6 represent our best approximations to 
this diffuse ISM cloud. Physically, the model distribu- 
tions represent cool, low density clouds surrounded by 
warm, lower density envelopes. Within the clouds, cold 
dense cloudlets are interspersed here on a regular grid. 
Cloudlet formation is beyond the scope of this paper but 
is likely a thermal, rather than gravitational condensa- 
tion. 

Model 5, shown in Figure|S](top), is of a Mach 10 shock 
over-running a cloud with a density distribution given by 
Equation^ The cloudlets have a x of 10 and a maximum 
extent of a c i ouc ii e t = 0.05. Furthermore, the inner radius 
of the cloud core is 0.6ao, which results in some of the 
cloudlets being outside of the cloud core. 

As seen in Figure El the blast wave shock is hardly 
slowed by it's interaction with the cloud, similar to Model 
4. However, the high density cloudlets do significantly 
alter the cloud shock structure. Compared to Model 2 
(Fig.|3J), the cloudlets appear to play a significant role in 
slowing the cloud shock. 

Model 6 differs from Model 5 in three ways: First, the 
X of the cloud is 8, rather than 10; secondly, the x of 
the cloudlets is increased to 15, and thirdly, the radius 
of the cloudlets is 0.03. Model 6 is shown in Figure [5] 
(bottom). Here one sees that the x °f the cloud is so 
low that it barely slows the shock. Moreover, and prob- 
ably more importantly, the spacing of the cloudlets is 
such that they do not feel the effects of their neighbors 
(i.e. Ar > 4.2a c z 011 <iZet & n d the diffracted shocks are not 
significantly curved). 
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Several features appear in the model simulations which 
are not observed. Prominent in all the models is the for- 
mation of a bow shock behind the cloud. A bow shock 
is not seen in Figure ^ simply because it is moving back 
into previously ionized material, so that there is no neu- 
tral population to excite. In models containing cloud 
clumps (Models 3-6), fingers and mushroom heads are 
prominent in the post shock flow. Three-dimensional 
models for shock cloud interactions show that many of 
these features are unstable and will not persist in three 
dimensions due to turbu lent effects in the post shock flow 
jStonfi fc NormanllT99^. 

6. DISCUSSION 

As seen in Figure^] the southwest cloud of the Cygnus 
Loop represents a fairly uncomplicated case for investi- 
gating many of the basic phenomena of a shock-cloud 
interaction. The low density of the cloud implies that 
the cloud shocks will be largely non-radiative in na- 
ture. Compared to other regions of the Cygnus Loop 
ijLevenson fe Grahamll200ll and references therein), the 
low density nature of this small cloud allows us to view 
internal cloud shocks. Furthermore, while the east-west 
extent of the cloud is not known, the location of the for- 
ward shock not interacting with the cloud suggests that 
this shock-cloud interaction is relatively young. Thus, 
this cloud presents a good test case to model the inter- 
action between a SNR shock and an ISM cloud. 

6.1. Comparisons to Other Shock Models 

There have been several previous studies on shock- 
cloud interactions. Perhaps the best current model for 
comparison is that of KMC94. Though there is not a 
one-to-one comparison between our Model 1 and their 
models due to the differing initial condition, many of 
their conclusions are observed in Model 1 such as the 
formation of K-H instabilities on the order of t cc (top, 
Panel c, Figure EJ), like that found in KMC94. 

On the other hand, there have been few published 
studies concerning the interaction between a shock and 
a cloud with a smoothly varying density. Our models 
seen in Models 2, 5, and 6 suggest that much of the in- 
stability growth observed in previous studies is related 
to the chosen geometry. ISM clouds are often modeled 
as cylinders or spheres with sharp, well defined bound- 
aries. Yet, the real boundary between the ISM and em- 
bedded diffuse clouds is likely to be less distinct. How- 
ever, models where the density varies over a certain dis- 
tance such as that described by a hyp erbolic tangent 
ijPoludnenko. Frank, fe Blackmanl [2002) can sometimes 
lead to spurious instability growth like that seen in the 
sharp edged cylindrical case. 

In regard t o the internal cloud density s tr ucture 
('cloudlets'), iPoludnenko. Frank, fe Blackmanl l)2002fl 
found that the principle parameter is the spacing be- 
tween the cloudlets. Their models suggested that 
there exists a critical separation between cloudlets of » 
4:.2a c i ou diet, and not surprisingly, the cloudlet spacing in 
Model 5 is about this value. Furthermore, a s pointed 
out bv IPoludnenko. Frank, fc Blackmanl l)2002ft . a larger 
Xcioudiet combined with a larger cloudlet spacing does 
not result in dynamics which are similar to the case of a 
lower Xcioudiet combined with a smaller cloudlet spacing. 
Instead, as evidenced by Model 6, the larger separation, 



regardless of Xcioudiet , results in what are essentially mul- 
tiple, independent interactions between the cloudlets and 
the cloud-shock. 

6.2. The Southwest Cloud 

While the southwest cloud represents a valuable lab- 
oratory for investigating the shock-cloud interaction, as 
evidenced in Figure Q it is still highly complex on small 
scales. Hence, the models presented here only approxi- 
mate its global properties. Based on Figure the cloud 
has a radius (N-S direction) of 1-2'. At the assumed dis- 
tance of 550 pc, this corresponds to 0.16-0.33 pc, or « 
0.5-1 x 10 18 cm. In Model 5, the fiducial radius of the 
cloud is 0.35. Using Model 5 as our potential model for 
the southwest cloud, the length scale of the model is thus 
1.5 x 10 15 cm cell -1 . 

The density of the ISM in this region is estimated to 
be ~ 0.1 - 0.3 cm~ 3 . The maximum x f° r Model 5 is 
10, but in reality the density profile is truncated at an 
inner radius = 0.21ao. At = 0.21ao, x ~ 4 f° r 
Model 5. This agrees with our lower density estimate of 
X ~ 4.5 from ram pressure arguments. Thus, the cloud 
particle density is w 0.4 - 1.2 cm -3 with a % = 10 for the 
individual clumps in Model 5 [n « 1.0 - 3.0 cm -3 ). The 
lower shock velocities seen in the cloud suggest cloudlet 
X's as high as 17, but the difference between a Gaussian 
profile with a peak density of 10 and one of 17 is minimal. 

Based on the size of the grid cell and the shock veloc- 
ity, the ambient shock traverses one cell in ~ 10 8 s ~ 3 
yr. The time difference in the proper motion analysis is 
about 10 years; that is, the ambient shock has traveled 
3-5 cells. In Figure El the top panels show the density 
at t = 2.2 - 7.3 x 10~ 2 . The simulation begins at t = 
0. and the shock first hits the cloud at t = 0.005. The 
radius of the cloud is sa 500 cells. Therefore, the ambient 
shock has been traveling for ~ 2000 yr when it reaches 
the cloud midpoint. 

Fr om the X-ray d e rived shock velocity of ~ 300 km 
s -1 . iPatnaude et all lj2002f) estimated the age of the in- 
teraction to be ~ 1200 years, so the modeled cloud size 
and the shock velocity appear reasonable. At the cur- 
rent epoch, the forward shock is 1' - 2' ahead of the 
cloud shock. From Figure the cloud shock lags behind 
the blast wave by 10% (bottom, Fig. Panel 6). This 
corresponds to a physical distance of 1.9 x 10 17 cm, or 
0'.5 at a distance of 550 pc. By Panel c of Model 5, the 
blast wave is 20% farther along than the cloud shock. 
Here, the morphology of Model 5 closely matches that of 
the southwest cloud (Fig.QJ. 

The observed internal cloud structure in the Ha im- 
age (Fig. is not that unlike the modeled shock cloud 
internal structure seen in Figure^ In general, the cloud- 
shocks seen in the Ha image are < 0f5 tip to tip. This 
scale is consistent with the approximate size of the inter- 
nal shocks seen in our Model 5 (Fig.|5j|. The cloud shocks 
have survived the 10 years between observations. The 
models, however, show that internal shocks are straight- 
ened out over a course of a few hundred time steps (~ 
200 yr). However, over the short time we are concerned 
with here, the shock structure of the cloud shock looks 
remarkably similar to that of the southwest cloud. 

7. CONCLUSIONS 
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A relatively isolated, low-density ISM cloud situated 
along the southwest limb of the Cygnus Loop provides a 
particularly clear view of the early stages of a SNR shock 
- ISM cloud interaction. The combination of multi-epoch 
observations and high resolution numerical modeling of 
this cloud has provided some new insights regarding how 
shocks overrun ISM clouds. The southwest cloud's isola- 
tion and low-density has also allowed us to view its in- 
ternal density structure and make inferences concerning 
the cloud's initial density distribution. 

Our specific findings are: 

1) Using multi-epoch Ha observations of a small, iso- 
lated ISM cloud in the southwest portion of the Cygnus 
Loop, we measured proper motions of Balmer-dominated 
shock filaments which wrap around the cloud, as well as 
the proper motion of several internal cloud shocks. The 
Balmer-dominated filaments have transverse velocities of 
~ 200-250 km s" 1 , while the shock filaments internal to 
the cloud have transverse velocities of 65 - 140 km s _1 . 

2) The shocked cloud's morphology does not show 
many of the dynamical instabilities predicted by previ- 
ous shock-cloud models. This suggests that there is no 
abrupt boundary or edge for diffuse ISM clouds. A sharp 
density rise between the cloud and the ISM would lead to 
steep velocity gradients at the shocked cloud - shocked 
ISM interface. These steep gradients would in turn lead 
to the onset of Kelvin-Hclmholtz instabilities, which are 
not observed. This conclusion contrasts with the shock- 
cloud interaction seen in the southeast of the Cygnus 
Loop, where the blast wave is thought to be interact- 
ing with a large, dense cloud, and instability growth is 
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clearly seen along the cloud-shock boundary. 

3) Our model hydrodynamic simulations suggest that 
ISM clouds are best modeled as a constant or smoothly 
varying core density embedded in lower density envelope 
which tapers to the surrounding ISM. Ram pressure equi- 
librium arguments suggest a cloud-ISM density contrast 
for this cloud of % = 5 - 17, with the lower corre- 
sponding to the diffuse regions of the cloud and the upper 
limit of 17 corresponding to the dense cloud clumps. 

4) A definite spacing of dense, small "cloudlets" inside 
the cloud is needed to generate the cloud's internal mor- 
phology as seen in the Ha image. As pointed out by 
IPoludnenko. Frank, fc Blackmanl l|2002(K clumps spaced 
too closely together interact with the shock as if they 
were one large clump, while those spaced too far apart 
behave as a set of individual clouds. Our models are 
consistent with the optimal spacing of d rr it ~ ^a c ioudiet 
ijPoludnenko. Frank. &: Blackmanl l2002|) . The observed 
internal cloud shock diffraction caused by these cloudlets 
is a short lived phenomena. As the cloud shock interacts 
with the cloudlets, the diffracted shocks re-order them- 
selves on a time scale of order a few cloudlet crossing 
times. 
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Simulations of Shock-Cloud Interactions 



TABLE 1 

Measured Proper. Motions and Estimated Shock Velocities 



Non-radiative Filaments Cloud Filaments 



Filament 




A 1 






Filament 










Region 


(mas) 


(mas yr — 1 ) 


(km s 




Region 


(mas) 


(mas yr ) 


(km 




- 1 ) 


Fl 


630 ± 55 


55 ± 5 


140 ± 


10 


CI 


620 ± 170 


55 ± 15 


140 


± 


30 


F2 


610 ± 55 


55 ± 4 


140 ± 


10 


C2 


325 ± 55 


30 ± 3 


80 


± 


10 


F3 


745 ± 55 


65 ± 3 


170 ± 


10 


C3 


575 ± 55 


50 ± 4 


130 


± 


10 


Fl 


665 ± 55 


60 ± 3 


155 ± 


10 


C4 


410 ± 55 


35 ± 2 


90 


± 


10 


F5 


810 ± 55 


70 ± 5 


180 ± 


10 


C5 


400 ± 55 


35 ± 5 


90 


± 


10 


F6 


910 ± 80 


80 ± 7 


210 ± 


15 


C6 


500 ± 55 


45 ± 5 


120 


± 


10 


F7 


1075 ± 110 


95 ± 9 


250 ± 


20 


C7 


340 ± 85 


30 ± 6 


80 


± 


15 


F8 


960 ± 80 


85 ± 6 


220 ± 


15 


C8 


550 ± 85 


50 ± 7 


130 


± 


15 


F9 


780 ± 55 


70 ± 2 


180 ± 


10 


C9 


380 ± 160 


35 ± 13 


90 


± 


30 


F10 


1025 ± 55 


90 ± 7 


235 ± 


10 


C10 


425 ± 130 


40 ± 11 


105 


± 


25 


Fll 


1010 ± 110 


90 ± 10 


235 ± 


20 


Cll 


280 ± 85 


25 ± 8 


65 


± 


15 


F12 


900 ± 140 


80 ± 12 


210 ± 


25 


C12 


295 ± 55 


25 ± 2 


65 


± 


10 


F13 


1110 ± 55 


100 ± 4 


260 ± 


10 


C13 


380 ± 55 


35 ± 4 


90 


± 


10 


F14 


1120 ± 85 


100 ± 6 


260 ± 


15 


C14 


410 ± 110 


40 ± 10 


105 


± 


20 


F15 


745 ± 55 


65 ± 4 


170 ± 


10 














F16 


625 ± 55 


55 ± 4 


145 ± 


10 














F17 


800 ± 55 


70 ± 3 


185 ± 


10 














F18 


715 ± 55 


60 ± 2 


155 ± 


10 














F19 


690 ± 55 


60 ± 4 


155 ± 


10 














F20 


720 ± 110 


60 ± 9 


155 ± 


20 














F21 


660 ± 55 


60 ± 4 


155 ± 


10 















a 1992.6 - 2003.7 

b Shock velocities assume a distance of 550 pc I Blair ct al. 2005) 
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TABLE 2 
Model Parameters 



Model 


tt() 


\ inn .i 


M 


^cloudlets 


Xcloudlets 


^cloudlets 


1 


0.30 


6 


20 








2 


0.25 


10 


10 








3 


0.30 


6 


20 


13 


10 


0.03 


4 


0.30 


3 


20 


13 


10 


0.03 


5 


0.35 


10 


10 


15 


10 


0.05 


6 


0.35 


8 


10 


15 


15 


0.03 



a Models 1, 3, and 4 have cylindrical density distributions. Models 
2, 5 and 6 have distributions corresponding to Equation HI 
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Fig. 2. — Left: One-dimensional shock profiles (top) and cross-correlation function( bottom) for region F14. The star appearing in the 
1992 data (top, extreme left) is the result of using a wider filter during the 1992 observations. Right: Cloud-shock profile (top) and 
cross-correlation function (bottom) for region C9. 




Fig. 3.— Density plots of Model 1 (top) and Model 2 (bottom). Model 1 is shown at t = 1.4, 2.6, 3.6, and 4.6 X 10~ 2 . Model 2 is shown 
at t = 3.7, 5,5, 6.9 and 8.2 X 1(T 2 . 
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Fig. 5. — Same as Figure^Jfor Model 5 {top) and Model 6 {bottom). Model 5 is shown at 2.3, 4.6, 6.6, and 7.5 X 10~ 2 . Model 6 is shown 
at 2.3, 4.3, 6.4, and 7.4 X 1CT 2 . 



